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We present a theoretical study of the director fields and energetics of nematic liquid crystal shells 
with two pairs of surface defects. The pairs of defects can undergo abrupt transitions between a 
configuration of maximum separation to at state in which the defects are confined to the thinnest 
hemisphere. We construct a phase diagram that maps out the stability and coexistence of these two 
configurations as a function of shell thickness and thickness inhomogeneity. Our results compare 
favorably with the experimentally observed transitions in nematic double emulsion droplets and 
explain their hysteretic character. 



INTRODUCTION 



Many systems in condensed matter physics and elas- 
ticity can be treated as two-dimensional, though only 
very few, like graphene [1^ and colloidal crystals at liquid- 
liquid interfaces [IHl], are truly monolayers. The theory 
of plates and membranes [5 , superfluid ^ and liquid 
crystal films [7 can all be neatly described by a reduc- 
tion of the number of spatial dimensions from three to 
two by assuming that the thickness is small compared to 
the other two dimensions and approximately constant. 
This reduction of dimensions usually simplifies the anal- 
ysis significantly, because the number of variables to solve 
for is reduced and in addition one can employ well de- 
veloped mathematical machinery such as complex anal- 
ysis that is well suited to tackle two-dimensional prob- 
lems. For instance, the use of conformal mappings has 
been applied succesfully in superfluid films [6 . Another 
example of the use of a conformal mapping was in the 
study of orientational order on a spherical surface [8 . It 
was found that in the ground state four defects of charge 
one-half reside on the vertices of a tetrahedron, inspiring 
the idea of self-assembly of liquid crystal coated parti- 
cles into a diamond structure [9 , which in turn triggered 
a tremendous research activity [IQH28]. Since the de- 
fects form very distinct regions on sphere, they can be 
functionalised chemically, an idea that has been realised 
for a similar system, namely metal nanoparticles coated 
with a monolayer of tilted molecules [29 . Experimental 
model systems of spherical nematics have revealed that 
the thickness plays a crucial role [17] . Besides the tetrava- 
lent ground state which contains charge one-half defects 
lines, bivalent defect configurations with charge-one de- 
fects have also been observed. The reason for this is the 
finite thickness of the liquid crystal shell. The charge one 
defect lines can lower their energies by escaping in the 
third dimension, thereby leaving two pairs of point de- 
fects at the interfaces, as an alternative to splitting into 
two charge one-half defects. Furthermore, the nematic 
liquid crystal shells generated by the encapsulation of a 
water droplet by a nematic droplet are strongly inhomo- 
geneous in thickness due to buoyancy. Similar to effects 
on the buckling and folding of inhomogeneous solid cap- 



sules [30 , the inhomogeneity has pronounced effects on 
the mechanics of the liquid crystal. One of the strik- 
ing manifestations of the inhomogeneity are very abrupt 
confinement and deconfinement transitions in shells with 
two pairs of surface defects. The investigation of this 
phenomenon, and more generally the theoretical study 
of the director fields and energetics of inhomogeneous 
bivalent nematic shells, is the main concern of this ar- 
ticle. Although the thickness truly makes this a three- 
dimensional problem, we are able to use two-dimensional 
techniques such as conformal mappings to find an Ansatz 
for the director field in spherical shells. This method is 
presented in detail in section [TT[ In sections III, we study 
the homogeneous shells as a function of thickness, tak- 
ing into account the elastic anisotropics. In section [TVl 
we construct a phase diagram for inhomogeneous shells 
that maps out the stability and coexistence of the con- 
fined and deconfined configurations as a function of shell 
thickness and thickness inhomogeneity. Our findings are 
shown to be in qualitative agreement with recent exper- 
imental studies. Finally, the effect of elastic anisotropy 
on the deconfinement transition is briefly discussed in the 
concluding section [V| 



II. DIRECTOR FIELDS IN BIVALENT 
NEMATIC SHELLS 

The experimental system under consideration is a ne- 
matic double emulsion droplet: a nematic liquid crystal 
droplet of radius R that encapsulates a smaller water 
droplet of radius a, as depicted in Figure [T] In this way, 
we can create spherical shells of nematic liquid crystal. 
We define a thickness h = R — a of the shell. Since in 
general the displacement of the inner water droplet out 
of the centre of the nematic droplet. A, is nonzero, h 
should be thought of as an average quantity. A surfac- 
tant or polymer is added to the inner and outer water 
phases for two reasons. First of all, it stabilises the dou- 
ble emulsion droplet, because it prevents the inner wa- 
ter droplet to coalesce with the continuous water phase. 
Secondly, it anchors the nematic molecules parallel to the 
interfaces. In modeling this experimental system we will 
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FIG. 1: Schematic of a nematic double emulsion droplet of 
radius R. The inner water droplet of radius a is displaced by 
an amount A along the vertical direction, thereby making the 
top of the shell thinner. 



employ elasticity theory for nematic liquid crystals, in 
which one constructs a Frank free energy functional as 
an expansion in spatial distortions of the local average 
orientation of the molecules, i.e. the unit director field, 
n(x), that respect the symmetries of the nematic liquid 
crystal 



F[n(x)] = i j dV (i^i (V-n)'+i^2(n-Vxn)' 

+ K3 (n X V X n)^^ - K24 J dS • (nV • n + n x V x n) 

(1) 

provided that we assume that these deformations are 
small on the molecular lengthscale. Here, Ki, Ks 
and K2A are elastic constants measuring the amount of 
splay, twist, bend and saddle-splay deformations respec- 
tively. In most of the work presented, we will work 
in the one-constant approximation, in which the splay, 
twist and bend elastic constant are taken to be equal: 
K = Ki = K2 = K^. Furthermore, we discard the sur- 
face term, effectively taking 2i^24 = K. For a typical 
droplet size of 50 iim the anchoring energy is much larger 
than this elastic energy. Therefore, we can take the pre- 
ferred tangential alignment of the nematic molecules at 
the interface as a constraint, thus forming a boundary 
condition complementing the free energy. Our approach 
to minimising the free energy with respect to the direc- 
tor field, will be to find a realistic Ansatz given certain 
locations of the defects. By varying these locations for 
different shell geometries we obtain the energy landscape 
as a function of defect positions, thickness and thickness 
inhomogeneity. The technique we employ to obtain the 
Ansatz is the method of comformal mappings. With the 
inverse stereographic projection we can find an Ansatz 
for a director field in a homogeneous shell (section II A ) . 
Then, by using the electrostatic analogy we ca n exp and 
the Ansatz to the inhomogeneous case (section II B). An 
additional numerical minimisation takes care of the es- 
cape of the disclination lines in the third dimension. 



A. The inverse sterographic projection and the 
Ansatz for the homogeneous shell 

The Ansatz for the director n of the homogeneous 
bipolar shell, with two straight disclination lines along 
the z-axis, simply reads 



n (x) = cos a ^ + sin a 



(2) 



where and cj) are the unit vectors corresponding to the 
zenith, ^, and azimuthal, (/), angles respectively. Note 
that a is the angle over which n = ^ is rotated at each 
point on the sphere with respect to a orthonormal ref- 
erence frame. Thus, the director fieldlines for a = 
and a = 7r/2 correspond to the meridians and circles 
of latitude. To find the Ansatz for any other locations 
of the the disclination lines, however, we perform an in- 
verse stereographic projection. A director field in the flat 
plane, minimising the free energy, is projected onto the 
concentric surfaces of spheres with radii, 7^, varying be- 
tween a and i?, i.e. a < IZ < that fill up the shell. 
Hereby, angles are preserved, i.e. this mapping is confor- 
mal. This director field contains two charge-one point de- 
fects, as we wish to eventually construct an Ansatz with 
two charge one line defects spanning the shell. Since the 
Euler-Lagrange equation, valid everywhere except at the 
defect cores, for the angular field, <l>, defined as the angle 
of the director with a cartesian reference frame {u^v)^ is 
Laplace's equation, it obeys the superposition principle. 
Therefore, the director field can be written as the sum of 
the director fields of two individual defects, and <l>2, 
positioned at ri and r2. Thus we write. 



<l> (r) = ^1 + <l>2 = a + CJi + UJ2, 



(3) 



where uji is the azimuthal angle in the reference frame 
that has as its origin, as is shown in Fig. |2^. Here, 
a is again a global constant. If we take = (iii,0) the 
defects lie on the ii-axis. We have 




FIG. 2: (a) Point defects (red dots) in the i^i;-plane located 
at {ui , 0) . $ at r is the sum of the single defect solutions 
$i = uJi. (b) Schematic of the resulting fieldlines. 



uji = arctan {v^ u — Ui) , 



(4) 



where the two-argument arctan (?/, x) is as the ordinary 
arctan (^), except that it takes into account in which 
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quadrant the point (x, y) is. The resulting director field, 
n = cos <l> u + sin ^ v. (5) 



with 



^ = a + ^ arctan {v^ u — Ui) 



(6) 



= 1,2 



is displayed in Fig. for a = 0. We rewrite this by 
substituting the following identities, 



u = coscj p — smcj cj, 
V = sin a; p + coscj a>. 



(7) 
(8) 



where uo and p are the azimuth angle and the radial dis- 
tance in the ii'U-plane, respectively. The result reads 



n = cos (<l> — cj) p + sin — uo) u). 



(9) 



This field is projected onto a sphere by means of an in- 
verse stereographic projection [8 , illustrated in Fig. 3] 
Each point on the sphere is represented by a point in the 



and 



= tan 1 - I sm ( 



(14) 



as the director field on the sphere, depicted in Fig. [4] 
In ref. [8^ it was shown that this field on the sphere min- 
imises the free energy in the one-constant approximation. 
It posseses two charge one defects at zenith angles 



9j = 2 arctan 



27^ 



(15) 



At the same time this expression is an Ansatz for a 






k27etan(^/2)»| 
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FIG. 4: Top view of the director field on the sphere, given by 
eq. (fTT). 



homogeneous shell with two straight disclination lines 
spanning the shell, provided we build it out of concentric 
spheres of radius IZ. The director lies along the spheres, 
including the special case that these spheres are the sur- 
face of the inner or outer droplet. Therefore, the tangen- 
tial boundary conditions are satisfied. 



FIG. 3: The inverse stereographic mapping given by eq. (10) 
of the the i/^v-plane onto the sphere with radius IZ. 



plane by the following relation: 
u -\- iv = 27^ tan 



(10) 



Then, noting that the planar polar unit vectors, (p, a>) 
are mapped to the spherical ones, {0, cj))^ we find 



n = cos ($-(/)) ^ + sin (<I> - (j)) 0. 
with <I> given by 

<I) = Q/ + ^ arctan (3^^, Xi) 



= 1,2 



with 



?ii = tan ( - 1 cos (j) — sgn {ui) tan ^ 



(11) 



(12) 



(13) 



B. An electrostatic analogy and the Ansatz for the 
inhomogeneous shell 

The concentric spheres that fill up the homogeneous 
shell are displaced if the shell is inhomogeneous. More- 
over, the disclination lines are no longer straight. To 
construct an Ansatz for the director in inhomogeneous 
shells we need to find equations for the displaced spheres 
and the defect lines. For this, we exploit an electrostatic 
analogy, namely, calculating the equipotential (solid in 
Fig. [5| and electric (dashed green) field lines of an in- 
finitely long charged line running parallel to a conducting 
plane (blue) at a distance d. By the method of images, 
solving this electrostatic problem is equivalent to solving 
for the equipotential and electric field lines of two oppo- 
sitely charged parallel running cylinders, or, equivalently, 
a 2D point charge and its mirror charge. These can be 
extracted from the complex potential [33[ [34] 



{w) — log 



id 



w — id 



(16) 
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FIG. 5: Equipotential (solid) and field (dashed green) lines of 
an infinitely long charged line, indicated with a cross, running 
along the y-direction, parallel to a conducting plane (blue) at 
a distance d. The two equipotential circles, drawn in bold, 
correspond to two non-concentric droplets whose centers, in- 
dicated by dots, are displaced by A. 



where rja and tjr are the vertical displacements from the 
origin of the inner and outer spheres. We take r] < 0, 
such that A > 0. This implies that the thinnest part 
of the shell is at the top, like in Figs. [5] and [l] The 
other spheres fill up the shell. Since the spheres are the 
surfaces of revolution of the circles around the z-axis, we 
obtain the equation for the spheres simply by addition of 



y to the left hand side eq. (19): 



(22) 



Two independently chosen electric field lines will serve as 
disclination lines. These lines run perpendicular to the 
equipotential lines, and thus perpendicular to the surface 
of the inner and outer droplet, as is demanded by the tan- 
gential boundary conditions. Similar to the calculation 
of the equipotential lines, one can obtain the equations 
for the electric field lines from eq. ([l8| (appendix[A|). We 
find that the electric field lines are also circles: 



(23) 



where w = x -\- iz is a complex number. Note that ip is 
a conformal transformation, just like the inverse stereo- 
graphic projection is, mapping a region bounded by two 
non-concentric circles into a vertical strip (see appendix 
[A[ ). Thus, the level curves of the real and imaginairy parts 
of {w) are the equipotential and electric field lines, re- 
spectively. These two families of orthogonal lines, to- 
gether forming what is known as a isothermic net, read 
mathematically 



{w)] = log 



w -\- id 



w — id 



= constant. 



(17) 



yj _|_ 

(w)] = arg ( ) = constant. (18) 

^ w — id ' 



We see that eq. 17 describes circles of Apollonius (see 
appendix |A| with inverse points ±id. We can rewrite 
them as 



+ (Z - T])'^ 



(19) 



with radius, TZ, and displacement with respect to the 
origin, rj: which are related by 



V 



(20) 



Rotating the circular equipotential lines around the z- 
axis creates the non-intersecting spheres. By choosing 
two of these spheres (drawn in bold in Fig. [5| as the sur- 
face of our inner and outer droplets in addition to a choice 
of we can fix the geometry. The relative displacement 
of the inner droplet with respect to its concentric postion, 
A, is given by 



(21) 



with radii, and displacements, e, now in the x- 
direction, which are related as follows: 



(24) 



Since only the circular arc that is inside the shell matters, 
we care about the points of intersection of the two defect 
lines with the spheres that fill up the shell. We would 
like to find the zenith angle on each sphere, Pi , that 
these points of intersection make. We assign a different 
character than 6^^, because pi does not have a constant 
value as it depends on the displacement (or radius) of 
the sphere. Let us therefore refine our definition of 9i as 
the zenith angle of the defect on the outer-most sphere. 
Now, the following geometrical relations hold: 



X = ±1Z sin Pi 

Z = T] -\-7l COS Pi 



(25) 
(26) 



Then, by substituting x and z in the eq. J23| and elim- 
inating d in favor of 1Z by applying eq. (|2Q|) we find an 
expression for as a function of Pi^ r] and 1Z (assuming 
sin ft 7^0): 



± 



1Z-\- T] COS Pi 

sin Pi 



± 



IZ -\- T] COS Oi 

sin Oi 



(27) 



where the last equality follows from the constantness of 
e^, as we are moving along the same circle. We find the 
solution for pi 



Pi = 2 arctan 



TL-ri 



(28) 



Not surprisingly, ft is increasing as the radius of the 
sphere is decreasing. If = or = tt, the disclina- 
tion lines are straight and ft = or ft = tt, respectively. 
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Next, we find 77 as a function of tiie spatial coordinates 
X, y and since it is tiie only variable, besides the para- 
metric dependence on d and on which j3i is depend- 
ing. Put differently, given some point in space, on which 
sphere is it? To answer this question we resort to eq. 
(22), yielding the following result: 



7?(x) 



r 



2z 



(29) 



We have now acquired all the necessary information to 
construct the Ansatz for the director field in an inhomo- 
geneous shell. We take the Ansatz for the director field 



in a homogeneous shell, eq. (11 ), and make the following 
replacements 



0^ ^ ft, 
Z ^ Z — T] 



(30) 
(31) 



The first substitution concerns the defect lines. The sec- 
ond accounts for the displacement of the spheres and im- 
plies the substitution 



disclination lines leaving a point defect at the inner and 
outer surface for each line (see Figs. ?? and|6|. We refine 
the mesh at these defects to obtain good accuracy on the 
rapidly changing director (fig. ??). 




FIG. 6: (a) A typical mesh used in the numerical minimisation 
is refined around the locations of the defects. (b,c) Director 
field of the shell for (b) ^ = 0.2 and (c) | = 0.77. The 
disclination lines escape in the third dimension leaving two 
pairs of boojums. Each pair is encircled in red in (b) The 
director field in (c) resembles a slightly distorted director field 
of a single nematic droplet. 



{Z-T]) 



x'^ -\- y'^ -\- {z — rf) 



(32) 



with j3 being the zenith angle on the displaced sphere. Fi- 
nally, together with eqs. (27)-(32) we obtain the Ansatz 



for the director in inhomogeneous shells with two charge 
one disclination lines: 



n = cos ($ - 0) /3 + sin (^ - (j)) 0. 
where ^ is now given by 

<|) = Q/ + ^ arctan (3^^, Xi) 



=1,2 



with 



and 



Xi = tan ( ^ I cos ( 



sgn (ui) tan 



yi 



tan I — I sm( 



(33) 



(34) 



(35) 



(36) 



The disclination lines can be put anywhere except for the 
south pole. In the case of a bipolar defect arrangement, 
i.e. 6i = and Oi = tt, we draw on each sphere the 
director given by eq. ([2|, with the substitution in eq. 
(32) and find an ansatz for the bipolar inhomogeneous 



shell that reads 

n (x) = /3 = cos (3 cos (/)x + cos /3 sin c 



sin/3z. (37) 



The Ansatz is then subjected to a numerical minimi- 
sation, employing the finite element method [35| suitable 
for non-trivial geometries, to ensure the escape of the 



III. ENERGETICS OF HOMOGENEOUS 
SHELLS 

First, we calculate the free energy for the Ansatz in 
([2|, in which the defects are located on opposite poles 
on the sphere, by integrating the free energy density over 
the spherical shell except for a cut-off region determined 
by < (9 < ^ and n - ^ < < ir with r the radial 
coordinate. The result reads 



F =27r {Ki cos^ a ^ Ks sin^ a) 



(38) 



R\og^-{R-h)\og ^^^^ ^) -2h] (39) 



+ 27r {Ks - 2K2a) h 
^^2^K U\og ^-^-{R-h) log - 2h 



(40) 
(41) 



Note that splay and bend deformations can be trans- 
formed into each other by tuning a, but the total energy 
is unchanged if Ki = i^s, as is shown graphically in Fig. 
[7^. Furthermore, note that the saddle-splay term is pro- 
portional to the thickness, in contrast to the twist and 
bend energy which both contain a logarithmic divergence. 
We therefore expect that, as a first approach, it is not 
so important in determining the defect locations. Since 
these defects repel each other, this director field is ex- 
pected to be the ground state. Upon placing the defects 
at a different separation from each other, eq. (11), leads 



to an increase in the elastic energv[8l-[TQ|. Moreover, the 
splay and bend cannot be efficiently transferred into one 
another by a global rotation (changing a), e.g. splay no 
longer vanishes for a = whereas it did for eq. ([2|. This 
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(a) 



■(b)- 



FIG. 7: The splay (solid blue), bend (dashed red), and their 
sum (dashed-dotted green) as a function of a when (a) O12 = vr 
and (b) 812 = O.Itt, in both cases | = 0.025. 



obstruction that makes the shell locally look like a slab in 
which the lines can escape. Rather, the point of view that 
a slight director distortion is induced by a single nematic 
droplet (resulting in an energy cost) is more appropri- 
ate in this regime. This cross-over is illustrated in Figs. 
?? and [6j So, for thin shells, it is energetically favor- 
able to decrease the distance between the boojums that 
form a pair. However, for thick shells, the opposite can 
be concluded: it is energetically favorable to increase the 
intrapair distance between boojums. Apart from the in- 
trapair interaction described above, there is an interpair 
interaction. This is repulsive in nature, thus irrespective 
of thickness, for homogeneous shells we always find the 
bipolar arrangement as the free energy minimum. 



is presented graphically in Fig. ^)p. Note that the direc- 
tor field minimising the free energy for Ki 7^ Ks is not 
equal to the Ansatz [H [HI HH HH [36]. Besides the elas- 
tic anisotropy the escape of the defect lines in the third 
dimension modifies the energetics. As a result, there are 
two pairs of boojums residing on the interfaces. We can 
effectively take the escape into account in our calculations 
of the energy by replacing the cut-off b by the thickness 
h and adding A.2t: Kh for each pair of boojums [IQl [37] . 
We obtain in the one-constant approximation 



F = 27rKLRlog 



2R 



(R-h) log 



2{R-h) 
h 



2.2h 



(42) 

In Fig. [8| we compare this analytical estimate with nu- 
merical results from our procedure outlined in the previ- 
ous section. We find a good agreement, in particular for 
small ^, as expected. In this regime, the free energy rises 




FIG. 8: Free energy of the bipolar shell as a function of thick- 
ness. The line is given by the analytical estimate in eq. (42) 



as the volume of the shell increases. For large ^, our re- 
sult deviates from eq. (42). Remarkably, as the thickness 



is increased, the free energy decreases after some critical 
value, I ^ 0.6. The size of the inner droplet, which is 2a, 
is no longer larger that the scale over which the escape 
happens, which is roughly h. As a result, when h becomes 
comparable to R^ the inner droplet no longer forms an 



IV. ENERGETICS OF INHOMOGENEOUS 
SHELLS 

A. Buoyancy versus elastic forces 

Before we study the effect of the thickness inhomogene- 
ity on the mechanics of the nematic liquid crystal, we first 
investigate its origin. In our experiments we observe that 
the inner water droplet is displaced along the vertical di- 
rection. This implies that gravity plays its part, but it 
does not necessarily mean that it is the density mismatch 
between the nematic and water that drives the motion of 
the inner droplet. Another possibility would be that the 
elastic forces push the droplet out of the center, while 
gravity only breaks the symmetry. To identify the ori- 
gin of the thickness ihhomogeneity, we will compare the 
magnitude of the elastic forces with Archimedes force. 
Therefore we map out the elastic energy as a function 
of the displacement ^ for several values of ^, as shown 
in Fig. [9] Our first observation is that the stability of 
this perturbation is a nontrivial function of the thick- 
ness. For ^ = 0.2,0.3,0.6 we observe that the energy 
decreases as a function of A. This is in agreement with a 
calculation in ref. [TT] done for ^ = 0.77. However, for a 
relative thin shell of ^ = 0.1 there is an elastic minimum 
for A = 0. Second, the magnitude of the elastic force 
is less than or of the order of fe ^ K ^ 10~^^7V. This 
is much smaller than the net force from buoyancy and 
the weight of the droplet = {pnem — Pw) gV with the 



volume of the water droplet V = ^ira"^. For a ^ 50 pm 
and a difference in density between 5CB and water of 
roughly 3 10 kg at room temperature [38 , we find 
fb^2 10~^^N. Therefore, we conclude that buoyancy is 
indeed responsible for displacing the inner water droplet 
from the center in our experiments. If one would try to 
match the density of the nematic to the water density, 
as was done in the experiments in ref. [11 , where the 
density difference was brought down to 2 kg , and 
fe will be of the same order, but only when the inner 
droplet is at the periphery. Also in the regime of small a 
these forces will become comparable. 
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FIG. 9: The elastic free energy of the Uquid crystal as a func- 
tion of the relative displacement of the inner droplet. In 
(a), (b) and(c),-^ = 0.6,0.3,0.2, respectively, the energy is 
minimised when the droplet is on the periphery of the larger 
droplet, resulting in an elastic force of the order of K push- 
ing the inner droplet outwards, (d) For a thinner shell with 
^ = 0.1, there is a restoring force on the inner droplet, driving 



it back to the center of the outer one. 




(b) 



(c) 



^2 



/tt 



^12, 



/tt 



a2 



I = 0.2 and f = 0.831 
= 0.8 and f = 0.946. 



FIG. 10: The free energy for a shell as a function of the angu- 
lar separation between the two defects on the outer surface in 
the confined configuration, when (a) 
(b) I = 0.7 and f = 0.939, (c) ^ 
The dashed line indicates the energy of the deconfined con- 
figuration. These graphs suggest a confined global minimum 
at ~ 0.17 in (a), a local confined minimum at ^ 0.1 and 
a global deconfined minimum in (b) and a global deconfined 
minimum in (c). 



C. Phase diagram 

We construct a phase diagram as a function of thick- 
ness and thickness inhomogeneity. We find that for a 
given thickness there is a deconfined minimum below a 
critical value of the relative displacement marked 



B. Confined and deconfined defect configurations 



green in Fig. 11, which is monotonously increasing with 
the thickness. The confined state is found to minimise the 



In the remainder of this article we will compare two 
distinct defect configurations. In one configuration the 
defects are at maximum angular separation from each 
other at opposite sites on the sphere. This we will refer to 
as the deconfined state. In the other case the defects are 
trapped or confined to the thinnest top part of the shell. 
The defects are located symmetrically at an angle Oi from 
the vertical axis so that their angular separation is simply 
0i2 = 26 i. The energy can be estimated to grow with the 
thickness of the shell where the defects are located. This 
is roughly the minimal thickness at the top of the shell, 
for which there is a simple geometrical relation hmin = 
h — A. From this one immediately sees that h and A 
take opposite roles. We thus expect the confined state to 
be energetically favorable over the deconfined state when 
the shell is sufficiently thin and inhomogeneous, i.e. low 
h and high A. This heuristic argument has led us to a 
systematic study of the energy landscape as a function of 
defect location. We classify three cases: I) the deconfined 

~ n) 



lOi 



state is the only energy minimum , see Fig. 
both the confined and deconfined state are minima, one 
of them is local and the other is global, see Fig. [TofD ; 
III) the confined state is the only energy minimum, see 

Fig. im. 



0.8 



0.6 



A 

h 



0.4 



0.2 



/ confined/ i 






deconfined 



0.8 



FIG. 11: Phase diagram of the confined (green), deconfined 
(purple) and coexistentance phase (blue) as a function of 
thickness of the shell, and thickness inhomogeneity ^ . The 
confinement (green), and deconfinement (purple), 
transition lines separate these phases. The red line represents 
the assumed experimental trajectory of constant hmin- 



energy above another critical value, -rp, marked in purple 
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in Fig. [TT] which is also larger for thicker shells. There- 
fore, as anticipated in the previous section, we find that 
the confined defect state minimises the elastic energy for 
thin and inhomogeneous shells, whereas the deconfined 
defect state minimises the energy for rather homogeneous 
and thick shells. Since these two critical values for ^ are 
different there exist two minima for ^ < ^ < We 
can thus divide the phase diagram in three regions: a 
deconfined minimum-only, confined minimum-only and 
coexisting region coloured purple, green and blue in Fig. 
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respectively. These phases are separated by lines 
marking where, as in a first-order phase transition, a lo- 
cal energy minimum is lost. We remark that the energy 
differences between the deconfined and weakly confined 
states for thin and homogeneous shells become too small 
to conclude with certainty that ^ goes to a finite value 



h 

and the deconfinement transition reaches 
tremely low h. 



= at ex- 




D. Comparison with experiment 



FIG. 12: Angular separation between the defects as a function 
of normalised thickness, in the experiment (black circles) 
and in theory (red squares). The open symbols depict the 
hysteresis. 



The nematic double emulsion droplets create inho- 
mogeneous shells, because buoyancy displaces the inner 
droplet upward from its concentric position along the 
gravitational direction. The short-range steric repulsion 
from the polymer polyvinyl alcohol (PVA), prevents the 
inner droplet from coalescing with the continuous phase. 
Therefore, we assume that the thinnest part of the shell, 
hmim is effectively constant. By osmosis the thickness 
inhomogeneity can be modified. We find 



= 133 nm, a= 125 



A 

1 



(i-«or 



(43) 



where u = ^ and uq is the value of u when the shell be- 
comes homogeneous, see appendix [B] This path through 
the phase diagram is indicated in red in Fig. 11 If we tra- 



verse this path in the direction of decreasing thickness we 
find that the angular separation between the defects, 6>i2, 
changes abruptly from tt to a value much smaller than 
that, as does the order parameter in a first-order phase 
transition. In the model this occurs in both theory (red 
squares in Fig. 12) at ii/i^o ~ 30 and in the experiment 



(black circles in Fig. 12) at u/uq ~ 20. The abrupt- 



ness of the confinement transition is marked by the the 
short timescale of only tens of seconds, compared to the 
hours over which the osmosis occurs, in which the pair of 
defects located at the thicker hemisphere moves toward 
the top of the shell (see Fig. 13). Upon decreasing the 



thickness and consequently the thickness inhomogeneity 
even further the defects spread and the angular separa- 
tion increases again. When the shell is approximately 
homogeneous (see Fig. 14a), the effect of confinement 
has weakened so much that the defects are aligned an- 
tipodally. The axis joining them can now point in any 
direction though, as shown by the two shells in Figs. 14b 




FIG. 13: Bright field images of a shell (top view) undergoing 
the confinement transtition. A pair of boojums moves from 
the lower hemisphere (top row) to the upper hemisphere (bot- 
tom row). The time span is tens of seconds. 2R = 133 /xm and 
2a — 125 iim. 



and 14c. In this case, the energy of the thin shell does 
not depend on the orientation of this axis, in contrast 
to what happens for thicker shells, whose boojums axis 
are aligned along the gravitational direction. This also 
confirms that the defect deconfinement transition in the 
phase diagram goes to A = for low h. Upon reversing 
the path through the phase diagram, i. e. traverse the red 



path in Fig. 11 in the direction of increasing thickness, we 
first find that the defects move toward each other gradu- 
ally. Upon increasing the thickness even further we find 
that 0i2 increases rapidly to its maximum possible value 
at ~ 50 and ^ ~ 80 in the model and experiment, 
respectively, as the mutual repulsion between the pairs of 
defects becomes too large. It is thus favorable to have one 
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FIG. 14: (a) Bright field image of two thin shells obtained 
after osmotically shrinking the corresponding thicker shells. 
The location of the defects are easily seen in cross-polarization 
by (b) focusing the shells at the top, to see the upper pair 
of boojums, and (c) at the bottom, to see the lower pair of 
boojums. The defects are a diameter away, with O12 = tt. 
Interestingly, for these thin shells, the direction of the axis 
joining the two pairs of boojum is not correlated with the 
gravitational z-axis. The dimensions of the upper shell are 
2R = 138.8 /im and 2a = 137.5 /xm, corresponding to u/uq = 
1.05. 



pair of boojums at the thickest part of the shell. Note 
that the thickness at which this deconfinement transi- 
tion occurs is thus larger than the thickness at which the 
confinement transition occurs. This hysteresis between 
the confinement and deconfinement transitions is due to 
phase coexistence. The green and purple curves in Fig. 
[TT| (corresponding to the confinement and deconfinement 
transition) intersect the red curve (assumed experimental 
path) at different points in the phase diagram. 



V. CONCLUSION 

In this study we have crossed from a two-dimensional 
description of a spherical nematic liquid crystal to a 
spherical bivalent shell with a finite thickness and possi- 
ble inhomogeneity. Irrespective of thickness, we always 
find an antipodal arrangement as the free energy min- 
imum in homogeneous shells of nematic liquid crystal. 
However, this scenario changes when the shell thickness 
is sufficiently inhomogeneous. The repulsion between the 
pairs of boojums competes with the minimisation of the 
distance between the defects within a pair. As a re- 
sult, the defects undergo a confinement transition to the 
thinnest part of the shell. Conversely, the defects con- 
fined in the thinnest hemisphere make a deconfinement 
transition that maximises their separation. The critical 
displacement of the inner droplet for which these transi- 
tions occur are in general not equal, i.e. there is hysteresis 
present. These transitions are also present in our experi- 
ment, where a water droplet encapsulates a nematic liq- 
uid crystal droplet to make a spherical nematic shell. We 
have showed that these shells are inhomogeneous due to 
the buoyancy that displaces the inner droplet along the 
gravitational direction. Additional to the confinement 
and deconfinement transitions, a continuous evolution is 
observed, when thin shells become less inhomogeneous. 
Though we found an excellent qualitative agreement be- 
tween theory and experiment for all these phenomena. 



an exact quantitative agreement is still lacking, possibly 
due to a lack of validity of the one-constant approxima- 
tion. It would be interesting to extend this study by 
investigating the role of elastic anisotropy on the defect 
transitions in nematic shells. Since it is more difficult to 
exchange splay and bend when the defects are confined, 
we expect that the region in the phase diagram occupied 
by the confined state will be smaller if elastic anisotropy 
is included. This would imply that the confinement and 
deconfinement transitions occur at smaller thickness if 
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Appendix A: Conformal mappings and the circles of 
Apollonius 

In this article, we make extensive use of the method of 
conformal transformations. To obtain the director field 
on the sphere minimising the free energy we use the in- 
verse steregraphic projection, in eq. (p^Ol). Another ex- 



ample of a conformal mapping is il) {w) in eq. (16) 



i/j {w) = log 



A 



B 



(Al) 



with —A = B = id and d real, to which we could asso- 
ciate electric potential to the real part, as this holomor- 
phic function must obey Laplace's equation. (Note that 
the analogy with electrostatics made in section |IIB| is 
to aid the explanation and not unique; we could have 
made an analogy with two-dimensional fluid flow just 
as well.) This maps an 'inhomogeneous annulus', i.e. 
the area bounded by two non-concentric circles, in the 
complex iL'-plane to a vertical strip in the complex tp- 
plane. Likewise, the non-concentric equipotential circles 
and electric field circles are mapped to vertical equipo- 
tential and horizontal electric fieldlines, as in a capacitor. 
It is the Mobius transformation 



r {w) 



w -\- id 
w — id 



(A2) 



that maps the inhomogeneous annulus to a homogeneous 
one, i.e. the region bounded by two concentric circles. 
Consequently, this annulus in the complex r plane is 
mapped to the vertical strip by the transformation 



= log r. 



(A3) 



The mapping is illustrated schematically in Fig. [Tsj 
The equipotential (solid black) and electric field (dashed 
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tp (w) 



rewrite eq. (|A5|) as 



arg {w — A) — arg {w — B) = 7, 



(All) 



In reference to Fig. [T6j consider two fixed points A and 
5 on a circle inscribing the triangle ABw^ where w is 
a third point somewhere on the circular arc AB^ then 
simple geometry tells u s the angle AwB, called 7, is con- 
stant. Now, eq. (All) simply follows from the fact that 



the sum of the angles of any (Eucledian) triangle should 
be TT. Alternatively, one can carry out the algebraic ma- 



FIG. 15: Schematic of the mapping of eq. (Al) 



green) lines in the w plane can thus be found by setting 
the real and imaginairy parts of ip {w) constant: 



{w)] = lo; 



w 



A 



w 



B 



arg 



w 



A 



B 



= constant, (A4) 
= constant. (A5) 



Eq. ( A4 ) describes the circles of Apollonius characterised 
by inverse points A and B. Instead of the more familiar 
specification of a circle as all the points that are a radius 
away from a center, eq. (A4) defines a circle as the locus 



of points for which the ratio of the distance to A and the 
distance to B is constant. It is straightforward to show 
that eq. (A4) indeed defines circles by rewriting it into 



x'^ -\- {z — df' 



(A6) 



where C is a constant. Some simple algebra now leads to 
the usual equation of a circle 



+ {z — rf) 



with radius, TZ, and displacement, rj: 



7^ 



1 + C 
1-C 

1 + C 



C 



d. 



(A7) 

(A8) 
(A9) 



A look at these eqs. yields the relation between the dis- 
placements and radii of the circles: 



(AlO) 



The electric field lines, which run perpendicular to the 
equipotential lines, are also circles. Since the argument 
of the product of two complex numbers is the sum of the 
arguments of the individual complex numbers, we can 













Vb 


arg [w — A) 


\ arg [w — B) 



FIG. 16: The electric fieldlines are circular arcs given by eq. 
(All), shown here for generic A and B. 



nipulations to find that 



(A12) 



with radii, and displacements in the x-direction, e, 
related in the following way: 



(A13) 



Appendix B: Experimental path through phase 
space 

In this appendix we derive the experimental trajectory 
through phase space, eq. (43), determined by a constant 
hmin- With a straightforward geometrical consideration 
we can relate hmin to A: 



^ = 1 

h h 



(Bl) 



As mentioned in the main text, the thickness is modified 
by a fiow of water with volume 5V that passes through 
the shell. We take SV > if the volume of the double 
emulsion droplet, Vtot^ is increased and 6V < if it is 
decreased. Upon writing SV = l^rv^, we obtain for the 
radius of the double emulsion droplet 



(B2) 



where R' is the radius of the double emulsion droplet 
before the flow of water. Since we assume that the volume 
of the shell, Vgheiu is conserved, the inner radius changes 
similarly 



a'3 + v\ 



(B3) 
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By invoking eqs. (B2) and (B3) we can write for 

. 3 



i-(i-^r 



(B4) 



We wish to find the path through phase space, that is, 
we want to write ^ as a function of ^. We find for the 



displacement (eq. (Bl)) 



A 



R' 



(B5) 



Now by using eqs. (B2) and (B3), note that 



h 



(B6) 



Substitution o f eg. ( |B4[ ) into eq. (B6), in turn substi- 
tuted into eq. (B5) yields 



A 

~h 



Uo 



{i-uY 



{l-uoT 



(B7) 



where 



with 



u = ^. We have chosen = I — uq^ 
Uo = 0.0153 the normahzed thickness at which the sheh 
becomes homogeneous. If we draw this trajectory in the 
phase diagram [TT] we observe that, as we decrease thick- 
ness, it crosses from the bipolar regime to the non-bipolar 
regime via the regime of coexistence. 
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